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We solve the 5 = 1/2 Kondo lattice model within the dynamical mean field 
theory. Detailed predictions are made for the dependence of the lattice Kondo 
resonance and the conduction electron spectral density on temperature and 
band filling, n c . Two low- energy scales are identified in the spectra: a renor- 
malized hybridization pseudogap scale T* , which correlates with the single-ion 
Kondo scale, and a lattice Kondo scale To <C T*, which acts as the Fermi- 
liquid coherence scale. The lattice Kondo resonance is split into a main 
branch, which is pinned at the Fermi level, and whose width is set by Tq, 
and an upper branch at lo ~ T* . The weight of the upper branch decreases 
rapidly away from n c = 1 and vanishes for n c < 0.7 (however, the pseudogap 
in the conduction electron spectral density persists for all n c ). On increasing 
temperature, we find that the lattice Kondo resonance vanishes on a temper- 
ature scale of order 10 To, the same scale over which the single-ion Kondo 
resonance vanishes in impurity model calculations. In contrast to impurity 
model calculations, however, we find that the position of the lattice Kondo 
resonance depends strongly on temperature. The results are used to make 
predictions on the temperature dependence of the low- energy photoemission 
spectrum of metallic heavy fermions and doped Kondo insulators. We com- 
pare our results for the photoemission spectra with available high-resolution 
measurements on YblnCu^ and YbAgCu^. The loss in intensity with increas- 
ing temperature, and the asymmetric lineshape of the low-energy spectra are 
well accounted for by our simplified S =1/2 Kondo lattice model. 

PACS numbers: 75.30.Mb,71.27.+a,75.20.Hr,79.60.-i 
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1. INTRODUCTION 



Heavy fermions exhibit a number of interesting ground states with 
anomalous low-temperature properties SHi'lfi!. These include paramagnetic 
heavy-fermion states down to the lowest temperatures (e.g. CeCu6, CeAbj, 
CeRu2Si2, YbAgCu4, and YbInCu4), magnetically ordered states with small 
ordering temperatures and reduced moments (e.g. URu2Si2, CeAl2), and 
unconventional superconducting states (e.g. CeCu2Si2, UBei3, UPta). This 
anomalous behavior is mostly attributed to the strong correlations between 
the electrons in the localized or weakly-itinerant / states of the rare-earth 
(or actinide) ions. The natural starting point for understanding this pecu- 
liar behavior is the periodic Anderson model or, when charge fluctuations 
are negligible, its low-energy counterpart, the Kondo lattice model. These 
models are fairly well understood in one dimension where extensive calcu- 
lations, based on exact diagonalizationLJ and density matrix renormaliza- 
tion groupS techniques are available. However, it is not clear how much of 
the insight gained is relevant to real- life three-dimensional heavy fermions. 
The Gutzwiller approximation, that can be applied in any dimension, has 
yielded some information on various correlated lattice modelsa, particularly 
for ground-state properties. The development of a non-pert uxbflt hre approx- 
imation scheme, the dynamical mean field theory (DMFT) li3E3l3l!Jli-3, has 
recently enabled a more comprehensive investigationof correlated lattice 
models, and particularly of their dynamic properties tU. In this paper, we 
apply this technique to the Kondo lattice model in order to investigate (i) 
the relevant low-energy scales of the model and how these appear in the 
single-particle and magnetic excitation spectra, and (ii) the temperature de- 
pendence of the /-electron spectra for the interpretation of the low-energy 
photoemission spectra (PES) ofnaramagnetic heavy-fermion compounds, of 
current experimental interestO'tZHls. 

The /-level photoemission and inverse photoemission spectra (IPES) of 
heavy fermions measure, respectively, the single-particle removal and addi- 
tion spectrum for / electrons. They provide information on the electronic 
structure of heavy fermions, including the energies of the relevant /-level con- 
figurations, their spin-orbit splittings, and the size of the screened Coulomb 
interaction between the / electrons. In addition to this information about 
the "high-energy" scales (the smallest of these being the spin-orbit splitting 
~ 1 eV), the /-level spectrum contains also information about low-energy 
scales associated with correlations and lattice coherence. Identifying these 
is one of the aims of the present study of the Kondo lattice model. 

We shall not deal, here, with the "high-energy" part of the spectra 
of heavy fermions, which is well accounted for by an Anderson impurity 
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model (AIM) approach llH^Hl. Instead, we shall focus our attention on 
the low-energy spectrum, which is the part most strongly influenced by lat- 
tice coherence effects, especially at low temperatures. These effects should 
emerge most clearly in the Kondo lattice model, which addresses specifically 
the low-energy scale, neglecting the complications associated with the high- 
energy atomic-like features. Although single-impurity model approaches also 
make predictions for the low-energy spectrum, there is no a priori reason to 
expect these to be quantitatively or even qualitatively correct. A lattice 
model calculation is therefore required to quantify the extent to which AIM 
calculations apply. 

In the absence of detailed lattice-model calculations, the experimen- 
tal results for the low-energy spectra of heavy fermions have, nevertheless, 
been discussed within the AIM approach. Two specific predictions of this ap- 
proach concern the magnitude of the temperature dependence of the near-^/p 
feature (Kondo resonance) in the spectra, and the scaling of the intensity 

predic- 
dispute 



of this feature with the Kondo scale Broad agreement with 
tions is claimed by some groupsll3'B3'Eac3E3 whereas otherstllcB 
the interpretation in terms of the AIM. The reader is referred to the re- 
views Ref. 17,1^, and to Ref. Pq^Gj JjC , 31 , for some recent aspects of this 
controversy. In the context of this discussion, these experiments address a 
very important issue, namely the temperature dependence of the low-energy 
spectra and the appropriate theoretical framework for interpreting them. We 
shall attempt to clarify this issue by calculating the temperature dependence 
of the lattice Kondo resonance and comparing this with the predictions of 
the impurity models and with experimental spectra. Experimental studies 
of this temperature dependence at temperatures both above and below the 
relevant Kondo scale are technically difficult, and so far have only been car- 
ried out on a few systems, for example on YbAgCu4c3, YbAl3c3, and more 
recently on YblnCuJill and URu2Si2 0. Most other studies have only com- 
pared spectra at two temperatures, typically room temperature and at one 
low temperature comparable to the Kondo scale: this is not sufficient for 
an unambiguous testing of model predictions with experiment, but can at 
most indicate agreement with predicted trends. Hence, it is important for a 
detailed comparison with theory for more experiments of the former type to 
be carried out. 

This paper is organized as follows: in Sect. we introduce the quantities 
of interest for the Kondo lattice model, and sketch the dynamical mean field 
method used; in Sect. |3] we discuss the main theoretical findings; in Sect. |4] 
we describe the results for the PES and IPES intensities; Sect. |T] contains 
the comparison of the theoretical spectra with actual PES data for YbInCu4 
and YbAgCu4; Sect. ^] summarizes the main findings and problems left open 
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by this paper; |Appendix A. provides details on the implementation of the 
DMFT to the Kondo lattice model and an expression for the self energy; 
Appendix B, defines the relevant spectral functions for this study. 



2. THE KONDO LATTICE MODEL AND THE DMFT 

METHOD 

2.1. The Kondo Lattice Model 

Consider the Kondo lattice model, given by the Hamiltonian 

k,cr j 

The first term describes non-interacting itinerant electrons with dispersion 
relation ek- The second term is the antiferromagnetic exchange interac- 
tion, of strength J, between the localized spins Sj at lattice sites j and 
the local conduction electron spin density Sj = Yl^, v c \ yPiivCjpi where 

1/9 

Cj >a = -/V sitcs ^ k Ck )(T exp(ik • Rj) destroys a local conduction electron or- 
bital at site j. This model is relevant for heavy fermions with negligible 
charge fluctuations on the / orbitals: it is appropriate for describing Ce 
and Yb compounds having a doublet crystal-field ground state for the 4/ 1 
(Ce) or 4/ 13 (Yb) atomic configuration, acting as a S = 1/2 pseudospin 
represented in the model by the Sj operator. 

For the unperturbed density of states of the lattice, po(u>), we take a 
semi-elliptic form 

PoH 6 k ) = ^^/W^ , (2) 

k 

corresponding to a Bethe lattice. The main features of the local dynamical 
quantities discussed below are primarily a result of the correlations and not of 
this particular choice of lattice, which is especially convenient in the DMFT. 



2.2. DMFT and Numerical Renormalization Group 

To treat the complicated correlation problem expressed by Eq. (|l|) we 
resort to a standard technique developed in the context of the infinite di- 
mensional limit of correlated lattice modelsS, the DMFtM In the DMFT 
approximation, the solution of (|lj) reduces to solving a quantum impurity 
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model embedded in an effective conduction electron medium, whose spec- 
tral density p(u) = ^ k 5(u> — ejj is determined self-consistently In 
the present case of the Kondo lattice model, the effective quantum impurity 
model is just the S = 1/2 Kondo modelE3 

Ti-imp = ^2 £k Ck, CT C k,(T + J s • s > ( 3 ) 
k,<r 

where S is the impurity spin and s is the spin density of the conduction 
electrons at the impurity site. The effective medium is determined from 
the self-consistency requirement that the local conduction-electron Green's 
function of the effective impurity model G c (lu,T) (which depends implicitly 
on ek) be identical to the local lattice Green's function Gl(u,T), i.e. 

G^u, T) = G Li<7 (uj, T)=J2 — ~~v~~? t\ ■ W 

^ oj + p - e k - 2J CT (u;, T) 

Here, T, a (uj,T) is the conduction-electron self energy, which is independent 
of k due to the local nature of the DMFT approximation0. In the following 
we drop spin indices for the Green's functions and the self energy since, 
in this paper, we shall only be interested in paramagnetic solutions of the 
Kondo lattice model. 

To solve the effective single-impurity problem (|3|) we use Wilson's nu- 
merical renormalization group (NRG)0 extended to finite-temperature spec- 
tra which is based on a logarithmic discretization that allows the 
low-energy scales (such as the Kondo resonance) to be accessed with high 
resolution. The self-consistency procedure involves: (i) guessing an initial 
form for the effective medium p(io) = X^k^( u; — ^k) 5 (h) diagonalizing (J3j) 
with the NRG and solving for its dynamics, (iii) re-calculating the effective 
medium by using Eq. (|j). Steps (ii)-(iii) are repeated until self-consistency 
is reached. Further details of the NRG and the DMFT procedure, which we 



refer to as DMFT(NRG), are given in Appendix A. For an application of 
the DMFT(NRG) to the finite-temperature Mott transition in the Hubbard 
model see Ref. |37|. 

In addition to the condition (Q) on the Green's functions, we also fix the 
number of electrons n c in the conduction band, by adjusting the chemical 
potential fi at each iteration, until n c converges to the desired value. This 
number is obtained in terms of the local conduction electron density of states 
p c (u,T), as 

/+oo 
dw f T (u)p c {u,T) , (5) 

3i a 
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where /t(^>) = 1/(1 + exp[(w — fj,)/T]) is the Fermi function (the Boltzmann 
constant and h are taken to be unity throughout the paper). The interacting 
conduction electron density of states, in turn, is defined by 

P c (lj, T) = -- Im G c (u, T) . (6) 

The chemical potential self-consistency is not implemented in impurity- 
model calculations, but it is clearly required for concentrated systems, such 
as heavy fermions. 

The quantity of main interest in this paper is the low-energy part of 
the /-electron spectral function of the Kondo lattice model A(uj,T). This 



is defined in Appendix B. It describes the Kondo resonance part of the /- 
electron spectral function in the vicinity of the Fermi level Ep. The PES 
(u < Ep) and IPES (u> > Ep) intensities are given in terms of A(u,T) 
by multiplying this by the occupation factors for electrons fr(u) and holes 
1 — /t(<^), respectively. Another useful quantity is defined in Appendix) 



B. : x'zzi^iT), the spectral density of magnetic excitations for the effective 
quantum impurity model. 



2.3. The Role of the Band-Filling n c 

We now discuss the significance of the other parameter in the Kondo 
lattice model: the band filling n c . The cj^ fermions represent renormal- 
ized Bloch states in the wide metallic bands of the crystal, with the correct 
symmetry for coupling to the localized / states. Even in a stoichiomet- 
ric compound, with integer number of electrons per unit cell, the effective 
occupancy of the c band need not be commensurate, since the conduction 
electrons are shared between the c band and other bands of different sym- 
metry, not included in the model. Thus, n c = 1 represents the special case 
of particle-hole symmetry in this effective band: it is appropriate for Kondo 
insulatorsfexamples of which include Ce3Bi4Pt3, SmBg, YbBi2, CeNiSn, 
and SmS E3), since p c , A and x'L show true gaps at T = 0. n c ^ 1 addresses 
naturally incommensurate filling of this effective band. In particular, n c ~ 1 
can be realized by domng the Kondo insulators, e.g. by substitution of Ce 
with La in Ce3Bi4Pt3E3 or of Yb with Lu in YbBi2c3. Intermediate values 
of n c ~ 0.5 describe metallic heavy fermions, which we study in detail in 
Sect. [fy. Finally, small values of n c should be relevant to low carrier-density 
heavy-fermion systems, e.g. Yb4As3 and YbBiPtR 

In addition to this direct significance of n c representing the actual num- 
ber of electrons coupled to the / states in different materials, one can also 
attach a second meaning to a variable n c : deviations of n c away from unity 
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Fig. 1. The spectral density A(lo,T = 0) of the lattice Kondo resonance 
(solid line) and of the conduction electrons p c (uj,T = 0) (dashed line), for 
(a) n c = 0.9, and (b) n c = 0.6. The detailed behavior near the Fermi level 
is shown in Fig. ||. 



introduce particle-hole asymmetry into the model. Such an asymmetry is 
generally present in the Anderson lattice model and results from an asym- 
metric position of the / level with respect to Ep. In the Kondo lattice 
model, the / level is projected out and real charge fluctuations are absent 
(n-f = 1). Although there is no / level in the latter, one can still partly 
simulate the effects of different /-level positions by changing instead the 
Fermi-level position Ep, i.e. by varying n c . Within this scheme, Ce heavy- 
fermion compounds correspond to rif + n c < 2, i.e. < n c < 1, and Yb 
compounds are the particle-hole counterparts of these with 1 < n c < 2. 
Here, we only compute spectra in the range < n c < 1, relevant for the 
spectroscopy of Ce systems: the complementary region 1 < n c < 2, relevant 
for Yb systems, is obtained by particle-hole symmetry exchanging lo — > — oj 
in the spectral functions. 



3. RESULTS 



All calculations are carried out for J/D = 0.3, sufficiently smaller than 
unity for the low-energy spectra to show universal scaling with the Kondo- 
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Fig. 2. The T = conduction electron spectral density p c (u, 0) (dotted 
line), the lattice Kondo resonance A(u, 0) (solid line), and the magnetic 
excitation spectrum x^' 2 (w,0) (dashed line), for decreasing band fillings: (a) 
n c = 0.9, (b) n c = 0.8, (c) n c = 0.7, and (d) n c = 0.4. The curves for \" are 
rescaled arbitrarily for the purpose of comparison with p c and A. Different 
frequency ranges are used in (a-d), as the lattice Kondo scale is decreasing 
with decreasing n c (see text). 



energy scale, but also sufficiently large to avoid unnecessarily large numbers 
of NRG iterations in the region of small band filling, where the Kondo scale 
becomes very small. All spectra are measured with frequency uj relative to 
the chemical potential p. 



3.1. T = Spectral Densities and Low-Energy Scales 

Figure |l| shows the zero-temperature spectra for p c and A on a wide en- 
ergy scale, of the order of the unperturbed bandwidth W = 2D, and for two 
fillings, n c = 0.9 and n c = 0.6. The large-scale effects of interaction on p c are 
clearly visible: band spreading, the opening of a renormalized hybridization 
pseudogap and consequently a reduction in height of the conduction elec- 
tron spectral density relative to the unperturbed spectral density (||). These 
global effects, on the scale of the whole band depend on the particular value 
of J and n c : we are interested instead in the low-energy behavior of the 
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spectral functions, which, instead, is universal, i.e. independent of J (but 
depend on n c ). 

Figure [2] zooms close to the Fermi level on the T = spectra, plotting 
them for different band fillings n c . Two different energy scales emerge clearly 
in these spectra. 

As soon as n c moves away from 1 (Fig. |2|a), the n c = 1 insulating gap 
turns into a pseudogap, for all spectra. We denote the characteristic width 
of this pseudogap in p c by T*. The pseudogap in p c lies just above the Fermi 
level for n c < 1 (and, by particle-hole symmetry, just below the Fermi level 
for n c > 1). It persists for all n c , but becomes increasingly asymmetrical and 
step-like as n c — > (see Fig. |2|b-d), making its precise value increasingly more 
difficult to specify accurately. Consequently, we have not made an exhaustive 
analysis of the dependence of T* on n c , but our calculations in the range 
0.2 < n c < 0.96 indicate that T* decreases slowly with decreasing n c and 
that it correlates approximately with the single-ion Kondo scale, being a 
factor of 2 — 3 larger than this, if we define the latter as the full width at half 
maximum (FWHM) of the single-ion Kondo resonance at T = 0. We find 
T*/D ~ 0.02 — 0.025 for n c > 0.6. The pseudogap in A manifests itself for 
n c close to 1. The ^4-spectrum represents the lattice Kondo resonance split 
by hybridization effects into a lower branch pinned at the Fermi level and an 
upper branch at uj ~ T*. The upper branch looses intensity at the expense 
of the main branch at the Fermi level with decreasing n c , and it eventually 
disappears for n c < 0.7. Our conclusion is that in the S = 1/2 Kondo 
lattice model, a splitting of the Kondo resonance by hybridization effects is 
only discernible close to half-filling (n c > 0.7), although in principle such 
a splitting is present for all n c (see Ref. |42"| , |43"1 ) . The pseudogap is present 
also in the spectrum of the magnetic excitations x'zz'- a peak at u; ~ T* for 
n c close to half-filling represents excitations across a spin pseudogap. On 
decreasing n c , the peak in x'L at « T* becomes a shoulder at n c = 0.8, 
and then it disappears for n c < 0.7, together with the upper branch of the 
lattice Kondo resonance in A. In contrast, for p c the pseudogap persists 
for all n c , so that T* remains a relevant scale for this quantity at all band 
fillings. 

The second energy scale apparent in the spectra of Fig. [2] is associated 
with the onset of Fermi-liquid coherence (as will become clearer with the 
discussion of the self energy below). It appears in Fig. ||| as the characteristic 
width of the T = lattice Kondo resonance in A (i.e. its main branch at 
the Fermi level), the lowest peak in the magnetic excitation spectrum x'zz 
and also as the energy range over which p c has a drop at the Fermi level. 
Clearly, there are several ways of defining this scale To. In this paper we 
define it as the position of the lowest-energy peak in the T = magnetic 
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co/T 

Fig. 3. Temperature dependence of the imaginary part of the self-energy 
Im[S(o;, T))] for n c = 0.8, showing that To is the Fermi liquid coherence 
scale. 



excitation spectrum x'z Z - The Kondo scale To, defined in this way, has the 
same dependence on n c and J as that defined in other ways, but may differ 
from these by factors of order unity (for example, the scale Tk defined as 
one half of the FWHM of the Kondo lattice resonance is 2=3 times larger 
than To). We have compared the dependence of To on n c with that of 
the single-ion Kondo scale, T; mp (defined in the same way as To, i.e., as 
the low-energy peak position in the local dynamical susceptibility, and using 
the same J and po as for the lattice model calculations). For n c > 0.4, we 
find that To < T; mp and To/Ti mp is consistent with the form n c exp(pra c ), 
with p ~ 5/2. found for the Anderson lattice model in the strong-correlation 
Kondo limitEl within the DMFT(NRG) approach. Burdin et a/El, using a 
slave-boson mean field approach, also find an increase of To/T; mp with n c . 

We find that, for fixed n c , both To and T* have the same J-dependent 
renormalization with decreasing J. To and T* are closest near half-filling 
and become increasingly differentiated with decreasing n c (see also Fig. ||||5| 
below). 
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Fig. 4. Temperature dependence of p c for n c = 0.8 in the region near the 
Fermi level and the pseudogap. The A-spectrum is shown for comparison. 
Temperatures and energies are shown in units of TqCj. 



3.2. T Dependence of the Self-Energy 

Figure || shows the temperature and frequency dependence of the imag- 
inary part of the conduction electron self-energy Im[E(o;, T)]. It vanishes at 
low temperatures at u = Ep on a temperature scale of order To. This is 
consistent with the interpretation of To as the Fermi-liquid coherence scale. 
The strong damping close to to ~ T* is associated with incoherent excita- 
tions due to the pseudogap in p c . It correlates with the position of the broad 
high-energy peak in the magnetic excitation spectrum of Fig. |2|b. 



3.3. T Dependence of p c and x'zz 

The role of the two energy scales T* and To in the spectra is further 
clarified in Fig. ||-[| which show the temperature dependence of p c for n c = 0.8 
(T*/T « 17) and n c = 0.4 (T*/T « 300). The temperature dependence 
of A is shown for comparison and will be discussed in detail below. We 
see that, as for the self-energy, the Fermi-liquid scale To also sets the scale 
for the temperature dependence of p c close to the Fermi level (just as for 
A). The feature of the p c spectrum directly affected on this temperature 
scale is the sharp drop at ^f- Instead, the region near the minimum of the 
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Fig. 5. Temperature dependence of p c for n c = 0.4 in the region near the 
Fermi level and the pseudogap. The A-spectrum is shown for comparison. 
The inset shows the pseudogap in p c on a larger scale. 



pseudogap in p c (and, whenever it exists, also in A and x'D starts to evolve 
at temperatures larger than To. It finally closes at higher temperatures 

For band fillings n c > 0.7, the pseudogap is also evident in x" z - Figure ^ 
shows that To sets the scale for the temperature dependence of a quasielas- 
tic peak, whereas excitations across a spin pseudogap (of order T*) give rise 
to a peak in x'lz at u; ~ T*. The temperature dependence of the latter 
is set by the scale T*. This type of response also characterizes the Kondo 
insulating behavior at n c = 1, as long as T > 0, so that the gap behaves ef- 
fectively as a pseudogap at finite temperature, as seems typical of correlated 
insulators^. The neutron scattering results for CeNiSnEJ and Ce3Bi4Pt3Ej 
show qualitatively the same temperature dependence. 



3.4. Discussion of the Two Energy Scales 

At high temperatures, approaching the scale T*, the spectra are similar 
to those of a Kondo impurity, namely the split Kondo resonance in A for 
n c > 0.7 reverts to a single- ion like Kondo resonance at the Fermi level whose 
height decreases with increasing temperature (see Fig. f?]a), x'zz nas a P e& k 



Temperature-Dependent Photoemission of Heavy Fermions 



15 



T/T, 



o' 



■9.78 



n=0.9 



T/T, 



o' 



6.46 



T/T, 



o' 



4.23 












20 



40 



60 



80 



co/T 

Fig. 6. Temperature dependence of x"z f° r n c = 0-9- This type of response 
is similar to the Kondo insulating behavior measured for CeNiSn e3 and 
Ce 3 Bi 4 Pt3§. 



at lo ~ T* whose height also decreases with increasing temperature, and p c 
is approximately independent of temperature. The physics appears like that 
of a collection of quasi-independent Kondo scatterers, since T* correlates 
with the single-ion Kondo scale describing the high-temperature properties 
of these Kondo scatterers. On decreasing the temperature, the effects of 
lattice coherence become quite pronounced for T rs T*\ the spectra develop 
a pseudogap resulting from the formation of hybridized bands between the 
/ and conduction states. Fermi-liquid coherence only sets in at much lower 
temperatures, specifically for T <^Tq. 

The formation of the paramagnetic heavy-fermion state therefore ap- 
pears to proceed in two stages, (i), a single-ion Kondo effect at T > T* 
in which the /-electrons are not strongly hybridized with the itinerant con- 
duction electron states, so that they play little role in screening, and, (ii), a 
low-temperature Kondo effect at T < T* in which the itinerant states have a 
strong /-electron character (via formation of renormalized hybridized bands) 
so that the latter also contribute to screening of the local moments. This 
second stage Kondo effect leads to the reduced lattice Kondo scale To de- 
scribing just the width of the main branch of the split Kondo resonance 
at the Fermi level. Two energy scales, a single-ion and lattice Fermi-liquid 
Kondo scale, have also been found within a recent slave-boson treatment of 
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co/T co/T 

Fig. 7. Temperature dependence of the lattice Kondo resonance for the 
same band fillings as in Fig. ^. Temperatures and energies are shown in 
units of To, where Tn is the lattice Kondo scale appropriate for each n c and 
defined in the text C3. 

the Kondo lattice modefl. 



3.5. T Dependence of the Kondo Resonance 

Figure [7| shows the temperature dependence of the lattice Kondo reso- 
nance for several band fillings. Compare the behavior of the main branch 
of the lattice Kondo resonance away from half filling with that calculated 
within the AIM 11 , and shown in Fig. ||[ The lattice Kondo resonance, like 
the impurity one, is pinned to just above the Fermi level for T <C To. How- 
ever, in contrast to the impurity case, the peak position E m (T) relative to 
the chemical potential shifts to below the Fermi level on a scale < To with in- 
creasing temperature. This effect is particularly evident slightly away from 
half-filling (e.g. n c = 0.9 in Fig. |?]a), and it becomes less pronounced as 
the band filling decreases toward zero, but is still significant as compared 
to the impurity- model predictions (Fig. ||). For the band fillings shown in 
Fig. 0, n c > 0.4, the lattice Kondo resonance vanishes on a scale of ap- 
proximately 10 To. This is the same scale over which the single-impurity 
Kondo resonance vanishes (see Fig. ||). We therefore do not see a signifi- 
cantly slower dependence (termed "protracted behavior"E2l) of the height of 



Temperature-Dependent Photoemission of Heavy Fermions 



30 



H 
3 



20 



10 



T/T. =11.41 

imp 

T/T =5.07 



T/T 
T/T, 
T/T h 
T/T. 



=2.25 
=1.01 
=0.45 
=0.20 



T/T imp =0.04 



1 ■ 1 1 1 1 1 1 1 ■ 1 

-20 -10 10 20 30 

Fig. 8. Temperature dependence of the impurity Kondo resonance calcu- 
lated from the AIM in the Kondo regime by using the The local 
level position used is = — 4A and the Coulomb interaction U = 4-7rA, with 
A = 0.1 -D the hybridization strength and D = 1 the half-bandwidth for a 
flat conduction band po(to) = 1/2D. The Kondo scale, T imp = 4.22 x 10~ 4 Z), 
is defined for the impurity model in the same way as for the lattice model, 
as the peak position in the T = magnetic excitation spectrum. 



the lattice Kondo resonance as compared to the impurity one. This con- 
tradicts the n c = 0.4 PES calculations for the periodic Anderson model of 
Ref 1511. Our calculations cannot rule out, however, that such "protracted" 
behavior might be present in other quantities, such as in the static suscep- 
tibility, or that such behavior will appear in spectral densities for n c < 1 
(the "exhaustion" limit!, an important conceptual limit that shall require 
clarification: we have not addressed this here, because we are mainly inter- 
ested in metallic heavy fermions with appreciable band fillings) . For a recent 
experimental evaluation of the "protracted screening" scenario see Ref. 53. 



4. T DEPENDENCE OF PES AND IPES 

The low-energy PES (/_) and IPES (i+) intensities are defined by 

I-(u,T) = f T (u>)A(u,T) (7) 
I+(u,T) = (l-f T (u))A(u,T). (8) 
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Fig. 9. Temperature dependence of the low-energy PES intensity (left) and 
IPES intensity (right) for the impurity Kondo resonance calculated from the 
AIM in the Kondo regime using the NRG K3. The same temperatures as in 
Fig. H| are used. 

For comparison with the Kondo lattice model predictions, we first outline 
the PES and IPES spectra from the AIM. 



4.1. T Dependence of PES and IPES for the AIM 

A typical T-dependent Kondo resonance from the AIM is shown in 
Fig. ||. The derived PES and IPES are shown in Fig. ^. It is often stated 
that the PES intensity in Ce compounds is very small compared to the IPES 
intensity because the Kondo resonance has most of its weight above the Fermi 
level. For the present N = 2 AIM we see that the Kondo resonance is, in 
fact, pinned very close to the Fermi level. Correspondingly, we observe nearly 
equal intensities in PES and IPES spectra in Fig. [9|. This does not contradict 
many PES and IPES experiments on Ce compounds, however, since most of 
these experiments are carried out at temperatures and resolutions far above 
the typical energies of the crystal field excitations or Kondo scale. As a 
result, they probe the full N = 6 degeneracy of the Ce 4/ 1 configuration 
and therefore cannot be described within an AIM with N = 2. Including 
the orbital degeneracy in the AIM is known to shift the position of the 
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Fig. 10. Temperature dependence of the low-energy PES intensity (left) 
and IPES intensity (right) for the lattice Kondo resonance for n c = 0.9 
(rather close to the Kondo insulating case n c = 1). The temperatures shown 
correspond to those in Fig. 0a. 



Kondo resonance to further above the Fermi levela, and therefore to increase 
the asymmetry between the PES and IPES intensities, as observed in many 
Ce (and Yb) heavy fermions. Nevertheless, in systems where the N = 2 
AIM approach is applicable, the spectra of Fig. |9] suggest that experiments 
at sufficiently high resolution and low enough temperatures (such that only 
the low-energy doublet is effectively playing a role) should give nearly equal 
intensities in PES and IPES for a system with a doublet ground state. 



4.2. T Dependence of PES and IPES for the Kondo Lattice 

Model 

The PES and IPES intensities for the Kondo lattice model are more 
complicated and show a range of behaviors, depending on band filling n c . 
Fig. [h]-12 shows the PES and IPES intensities for the four band fillings 



n c = 0.9, n c = 0.8, n c = 0.7 and n c = 0.4 of Fig. g. 

For Kondo insulators (n c = 1), there is experimental evidence on Ce3Bi4Pt3t 
and CeRhAsEl indicating the development of a gap in the PES on decreas- 
ing temperature. This is consistent with our predictions for the j4-spectrum 




Fig. 11. Temperature dependence of the low-energy PES intensity (left) 
and IPES intensity (right) for the lattice Kondo resonance for n c = 0.8. The 
temperatures are those of Fig. 0b. 



at 7i c = 1, to be published elsewhere^ and also with calculations on the 
symmetric Anderson lattice model e3>e1 

The results for n c close to 1 (Fig. 10) are relevant to doped Kondo 
insulators, e.g. to hole-doped Ce systems, such as Ce3_ x Lai_ x Bi4Pt3l2^l. As 
noted earlier (see Fig. ^), the Kondo resonance is highly asymmetric for 
n c < 1 with most of its weight lying below the Fermi level. Consequently, 
the bulk of the intensity in Fig. 10 lies in the PES part of the spectrum. In 
addition the strong shift of the Kondo resonance from above to below the 
Fermi level on increasing temperature for the case n c < 1, noted in Sect. [3] , 
further diminishes the near-Ep intensity in IPES as compared to that in the 
PES. 

Another interesting feature of our calculations is that the pseudogap can 
be seen as a dip in the IPES intensity at to ~ T*, which becomes increasingly 
well developed with decreasing temperature. Our predictions for hole-doped 
Ce systems apply equally well, by particle-hole symmetry, to electron-doped 
Yb systems. Results for the electron-doped Kondo insulator Ybi_ x Lu x Bi2^3 
indeed indicate an increase in the photoemission intensity at the Fermi level 
upon doping YbE>i2 with Lu. 

On further decreasing the band filling (Fig. [H|-[l2|), the IPES intensity 
increases with respect to the PES intensity. This results from the more 
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Fig. 12. Temperature dependence of the low-energy PES intensity (left) 
and IPES intensity (right), for the lattice Kondo resonance for n c = 0.4. 
The temperatures are those of Fig. [7|d. 

symmetric lattice Kondo resonance together with the tendency of this res- 
onance to move further above the Fermi level for smaller n c . In addition, 
the pseudogap in the IPES intensity fills in below n c ~ 0.7. For n c = 0.4, 
far enough from the Kondo insulating state to be representative of metallic 
heavy fermion systems, we find a behavior rather similar to the AIM spectra, 
i.e. the IPES intensity is slightly larger than the PES intensity at low tem- 
peratures, T <C To. However, a small increase of the temperature to T > To 
makes these intensities nearly equal (e.g. Fig. [l2] ). Inclusion of the orbital 
degeneracy of the /-level should lead to a displacement of the lattice Kondo 
resonance further above the Fermi level, as occurs in impurity models, and 
consequently to larger IPES intensities for Ce systems than PES intensities. 



5. COMPARISON WITH EXPERIMENTS 

Ideally, our predictions for the S = 1/2 Kondo lattice model are suitable 
for heavy-fermion systems with a local doublet ground state characterized by 
a separation Acf from the first crystal-field excitation much larger than the 
Kondo scale To. An example of such a system would be CeCu2Si2, which has 
Tq ~ 6 K and Acf ~ 30 meV Ea. Photoemission spectra have been measured 
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for this compound^ at 6 meV resolution (FWHM). The absolute resolution 
is high, but the relative resolution for probing the Kondo resonance is still 
quite low due to the small Kondo scale. Moreover, most measurements were 
taken at temperatures T > lOTo, with only one temperature (T = 15 K) 
relevant for testing predictionsabout the Kondo resonance. In this respect, 
the experiments on YbAgCuJ^I and YbInCu4Ell are more promising, since 
measurements were taken at T both above and below the relevant Kondo 
scales. Unfortunately, for these systems, neutron scattering studies £9'u indi- 
cate that the crystal-field excitation energy is comparable to the Kondo scale 
(Acf ~ To). In the temperature range of the experiments T ~ 10 — 200 K, 
thermal population of the crystal-field states will be important and would re- 
quire further refinements of our 5 = 1/2 model, which, however, are beyond 
the scope of the present work. With this proviso, we compare our results, 
which include lattice coherence and correlation effects, with measurements 
of the temperature dependence of the low-energy /-spectrum for YblnCuJ^ll 
and YbAgCu 4 ll. 

YbInCu4 exhibits an isostructural phase transition, similar to that of 
Ce, with a 0.5% volume change and a 0.1 change in /-level occupancy at 
a temperature T v = 42 K. Both phases show the characteristic features of 
paramagnetic heavy-fermion behavior. The high-temperature phase has a 
Kondo scale To = 30 K, whereas the low-temperature phase has a much 
larger Kondo scale To = 400 K 153. Temperature-dependent photoemission 
measurements were carried out in both phases^ (T =12-40 K for the low- 
temperature phase, and T =50-150 K for the high-temperature phase). 

YbAgCu4 showsa single paramagnetic heavy-fermion phase down to the 
lowest temperatures^. Neutron-scattering experiments indicate a quartet 
crystal-field ground state with a Kondo scale To = 100 K 0, consistent 
with thermodynamic measurements!^. Measurements for this compound 
were carried out in the temperature range 20-260 K and show a strong 
temperature dependence of the neav-Ep feature in qualitative agreement 
with expectations from Kondo phenomenologyil. 

The experiments on YbInCu4 used single crystals in contrast to the 
polycrystals used for YbAgCu4. In both experiments, one observes a near- 
Ep feature (due to excitations from the f 7 ^ 2 ground state of the Yb ions), 
interpreted as the Kondo resonance, and an / 5 / 2 spin-orbit side band at 
approximately 1 eV below Ep. The bulk and surface contributions to these 
excitations can be separated in the experiments, allowing bulk effects to be 
investigated. The / 5//2 spin-orbit side-band is not included in our model, 
but its absence will not affect the overall predictions for the / 7//2 excitation 
in a range of <C 1 eV around the Ep spectrum. 

The experimental resolution (AE = 25 meV FWHM) is included by 
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Fig. 13. Comparison of the calculated PES intensity for the lattice Kondo 
resonance with the experimentally measured intensity " for the J 7 / 2 near- 
er feature in YbInCu 4 at T = 12K < T = 400 K. The calculated curve 
includes the AE = 25 meV Gaussian instrumental resolution. 

convoluting the calculated intensities (described in Sect. ^]) with the appro- 
priate Gaussian broadening. Since the resolution is comparable to or below 
typical Deby e ener gies, some effect of phonon broadening on the spectra may 
be observablecSEll. Specifically, for both YbInCu4 and YbAgCu4, we take a 
Gaussian phonon broadening term of the form AE p h = aT reaching 50 meV 
at 150 K c3: the total width used to convolute the calculated intensities is 
then given by 

AE tot = {AE 2 ph + AE 2 ) 1/2 • (9) 

For all the comparisons, we use an intermediate value n c = 0.4. This 
choice is motivated by the requirement of being in a regime representative 
of metallic heavy fermions, characterized by an asymmetric position of the 
localized / states. 



5.1. Low- Temperature Phase of YbInCu4 



In Fig. 13 we compare the calculated and experimental intensities for 
YbInCu4 at the lowest temperature T = 12 K at which measurements were 
taken in the low-temperature phase (the curves are normalized so that the 
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peak heights coincide). Since the Kondo scale To = 400 K is much larger 
than the temperature, this measurement essentially probes the T = lattice 
Kondo resonance below the Fermi level, convoluted with a Gaussian of width 
AE = 25 meV. All other measurements in this low-temperature phase are 
for temperatures T < T v = 42 K <C Tq = 400 K and can also be regarded 
as essentially zero-temperature measurements. At such low temperatures, 
the Kondo resonance is fully developed, so no temperature dependence is 
expected, and indeed almost none is seen Ell. 

Figure 13 shows reasonable agreement between the calculations and 
the experimental results. In particular, the theory reproduces the correct 
width of the near-Ep feature, as well as the slow decay of the intensity at 
negative energies. We attribute the latter to the non-Lorenzian decay of 
the Kondo resonance at high energies similar to that found for the impurity 
Kondo resonance^ (whose lineshape for \uj\ ^> Tq is of the Doniach-Sunjic 
type). The position of the peak in the intensity at T <C Tq is related to the 
position E m of the Kondo resonance in A(u,T). From Fig. | we estimate 
that E m (T — > 0)/Tq —0.5 for the present calculations at n c = 0.4. The 
experiments would indicate a value closer to —1, i.e. larger in magnitude than 
the S = 1/2 Kondo lattice model result. The predictions of the AIM in the 
Kondo regime (with occupancy ~ 1) would be in far worse agreement, since 
from Fig. || we see that the Kondo resonance is pinned even more closely 
to the Fermi level than in the Kondo lattice calculations. It is important 
to note, however, that there is no fundamental relationship between E m 
and the Kondo scale. E m simply reflects a particle-hole asymmetry in the 
model, but on physical grounds one expects it to lie within approximately 
To of the chemical potential. Its precise value, both for impurity and lattice 
models, depends on specific details influencing the particle-hole asymmetry, 
such as the ground-state degeneracy of the / state, the presence of low- 
lying crystal-field states, and the amount of mixed valence as reflected in 
the /-level occupancy (fixed to nj = 1 in the Kondo lattice model)E3. As 
our model describes a local ground-state doublet with nj = 1, it would be 
surprising if it could capture the precise value of E m /TQ for YbInCu4, which 
would depend on such details. Hence the position of the peak in the intensity 
at T <C Tq is also not captured precisely. It is likely that a Kondo lattice 
model, including the /-orbital degeneracy, or the corresponding Anderson 
lattice model, will show complete agreement between theory and experiment. 
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Fig. 14. Comparison of the temperature dependence of the calculated PES 
intensity for the lattice Kondo resonance (lines) with the experimentally 
measured intensity (symbols) for YbInCu4 c5. The measurements were for 
the high-temperature phase, with a Kondo scale of To = 30 K. We com- 
pare only with the / 7//2 near-Ep feature (the _f 5//2 spin-orbit side band at 
— 1.375eV in the experiments is not included in the present model). The 
calculated curves are Gaussian convoluted to account for the instrumental 
and phonon broadening AE to t of Eq. (|9|) . 



5.2. High- Temperature Phase of YbInCu4 

We turn now to the high-temperature phase of YbInCu4, for which 
To = 30 K and measurements have been taken up to 150 K= 5 To El 
This range of temperatures is interesting, because the results of the previous 
section predict a strong temperature dependence in the near-Ep spectra in 
the range < T < 10 To due to the Kondo effect. A comparison of the 
measured intensities with the Kondo lattice model predictions is shown in 



Fig. 14. As in Fig. |13|, the resolution and Kondo scales are taken from 
experiment and no attempt was made to obtain a best fit to the spectra by 
adjusting To or AE. A phonon broadening reaching 50 meV at 150 K is 
included, as described above. This improves agreement for the two highest 
temperatures, but phonon broadening is not the main source of temperature 
dependence in the spectrum, which is instead the temperature dependence 
of the Kondo resonance in the spectral function A{uo,T). The experimental 
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Fig. 15. Temperature dependence of the low-energy PES intensity of 
YbAgCu4 c2l compared to Kondo lattice calculations. Phonon broadening 
and instrumental resolution are accounted for as for YbInCu4 in Fig. 14 



spectra are consistent with this interpretation, because they show a loss 
in intensity of the near-E'p feature as the temperature is increased from 
50 K to 150 K: this agrees with the loss in spectral weight of the Kondo 
resonance with increasing temperature, and cannot be attributed to phonon 
broadening, which conserves the total spectral weight. 

The lineshape is well reproduced by the calculations, as is the tem- 
perature dependence of the intensity. The curves are normalized so that 
the lowest-temperature peak heights in theory and experiment agree. This 
normalization gives agreement between theory and experiment for the tails 
and peak heights at all temperatures. Note that, at finite temperatures, 
the position of the intensity peak no longer reflects the Kondo resonance 
position (and therefore also not the Kondo scale). Thermal broadening of 
the spectral function A(uj,T) in J_, will shift the intensity peak to lower 
energies with increasing temperature, as observed in the calculations. The 
remaining discrepancy is essentially due to the oversimplified description of 
the /-states, resulting in too small a value for E m (T), as discussed above for 
the low temperature case. 
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5.3. YbAgCu 4 



Figure 15 shows a comparison of the present calculations to early ex- 
perimental results of Weibel et al. il for the temperature dependence of the 
near-£^F feature in YbAgCu4. As for YbInCu4, we take the experimental 
resolution and Kondo scale without any adjustment. The agreement is less 
good than for YbInCu4, but still consistent with the interpretation of the 
temperature dependence of the spectrum as arising mainly from the Kondo 
effect. The discrepancy in the peak position has the same origin as that 
described above for YbInCu4, and it could be overcome by a more realistic 
modeling of the / states of this compound. 

For this compound, we normalize the curves to the tails in the intensity 
at negative energies. That part of the spectrum has much weaker tempera- 
ture dependence, it reflects the non-Lorenzian tail of the Kondo resonance 
and it is the safest to use in case there are extrinsic (non-Kondo) temper- 
ature dependences close to Ep. In contrast to the case of YbInCu4, we do 
not obtain agreement for both the tails and the peak height of the intensity 
as a function of temperature. We conclude that these early data, although 
indicating a Kondo effect in the temperature dependence of the spectrum 
close to the Fermi level, are not fully consistent with our model. To investi- 
gate the origin of the discrepancy between theory and experiment, it would 
be interesting to repeat such experiments, particularly on single crystals as 
for YbInCu 4 0. 



6. CONCLUSIONS 

In summary, we have studied the temperature dependence of the single- 
particle spectral densities of the Kondo lattice model and their evolution 
with band filling within dynamical mean field theory. We found that the 
spectra exhibit two energy scales: a renormalized hybridization pseudogap 
scale T*, which evolves from the true hybridization gap of the Kondo insu- 
lating state at n c = 1 on moving away from half-filling, and a typically much 
smaller Fermi- liquid coherence scale To , characterizing the formation of the 
main branch of the split lattice Kondo resonance at the Fermi level. The 
weight of the upper branch of this lattice Kondo resonance at to ~ T* be- 
comes negligible for n c < 0.7. The scale T* was found to correlate with the 
single-ion Kondo scale. The high-temperature behavior of the local spectral 
densities is similar to that expected for a collection of independent Kondo im- 
purities. We found that the lattice Kondo resonance develops over the same 
temperature range, < T < 10To < T* , as in the single-impurity models, 
at least for the range of band fillings studied in this paper: 0.4 < n c < 0.9. 
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In this respect, the Kondo effect in this low-temperature phase, T < T* , ap- 
pears similar to the single-ion Kondo effect. However, the electrons involved 
have a strong / character arising from the renormalized bands formed below 
T ps T*. We interpret the above findings as indicating a two-stage Kondo 
effect for the formation of the paramagnetic heavy-fermion state: a single- 
ion Kondo effect at T 3> T* , and a low-temperature Kondo effect at T < T* 
resulting from the modified electronic structure below T*, and giving rise to 
the new scale To. Lattice coherence, associated with pseudogap formation 
on a scale T*, sets in at temperatures well above the Fermi liquid coherence 
scale To. 

The results for the temperature dependence of the lattice Kondo res- 
onance were used to calculate the low-energy part of the /-level PES and 
IPES intensities of heavy-fermion systems. We made detailed predictions 
for the spectra of doped Kondo insulators (n c < 1) and metallic heavy 
fermions. For the former, we found that most of the intensity of the Kondo 
resonance lies in the PES, and not in the IPES. This applies to hole-doped 
Ce or electron-doped Yb Kondo insulators. The predicted pseudogap in the 
/-electron spectrum is best seen in the PES of electron-doped Yb Kondo 
insulators. For hole-doped Ce systems, this would lie just above the Fermi 
level in the IPES spectrum. On moving away from the Kondo insulating 
toward the metallic heavy-fermion case, we found that the main intensity of 
the Kondo resonance shifts to the IPES spectrum and the PES and IPES 
intensities appear qualitatively similar to impurity-model calculations. 

We made a parameter-free comparison of our results with high-resolution 
PES data for YbInCu4 and YbAgCu4, taking the Kondo scales and resolu- 
tion from experiment. Good agreement was found for the lineshape and 
temperature dependence of YbInCu4, although we argued that a more real- 
istic treatment of the /-states (including, for example, orbital degeneracy) 
would be required to obtain the precise peak position in the PES intensity 
of the near-^p feature. In the experimental spectra we found evidence of 
some phonon broadening at higher temperatures, but the main temperature 
dependence of the near-Ep spectrum is due to the Kondo effect. 

The present study has focused attention on the paramagnetic metallic 
state of the 5 = 1/2 Kondo lattice model and has revealed some of the 
physics of this state, such as the relevant energy scales involved and how 
these appear in the PES and IPES spectra of heavy fermions. It would be 
interesting to identify these scales in other properties, such as in the optical 
conductivity, transport properties and in the thermodynamics. Work on this 
is in progress and will be reported elsewhere. 



Temperature-Dependent Photoemission of Heavy Fermions 

ACKNOWLEDGMENTS 

Peter Wolfle has played a key role in encouraging work by one of the 
authors (TAC) on correlated electron systems, and it is a great pleasure to 
dedicate this paper to him on the occasion of his 60th birthday. We thank 
the Newton Institute - Cambridge and the ILL - Grenoble for support during 
the development of this project. NM also acknowledges support from the 
SISSA-Trieste and Italian INFM. We wish to thank Y. Baer, N. Brookes, D. 
Moore and S. J. Oh for useful discussions. 



Appendix A. NRG AND DMFT DETAILS 

We provide, here, some of the details of the NRG and the DMFT applied 
to the Kondo lattice model. 

Starting from Eq.(|3|) we follow the standard procedure for casting the ef- 
fective impurity model into a linear chain form suitable for an NRG calculations. 
This first involves a logarithmic discretization of the conduction band about 
the chemical potential n(T) (which for the purposes of this Appendix is taken 
as the zero of energy). The discretized energies are ± D± A~ n , n = 0, 1, . . ., 
where ±D± are the upper (+) and lower (— ) band edges, respectively, and 
A > 1 is the discretization parameter, set to 1.5 in the calculations. A 
unitary transformation of these states then yields a new tridiagonal set 
b n ,a, n = 0, 1, . . ., such that 6o j0 - = co )CT - The Hamiltonian takes the following 
form: 



oo oo 

+ ^ AnOl.A+l,* + H - c -)- ( A -!) 



n=0,u n=0,u 



We use this form for the actual calculations, whichinvolve iteratively diag- 
onalizing Eq. (A.l) in the standard NRG manner a, i.e. adding the orbitals 



n = 0, 1, . . . successively and retaining only the lowest energy eigenstates at 
each stage. We retained 462 non-degenerate states at each NRG iteration 
and we checked that the spectra did not change significantly by further in- 
creasing the number of kept states. The spectra are calculated at the natural 
NRG temperature scales, T\r ~ AaIEI. 

This NRG technique is the basic ingredient in the solution of the Kondo 
impurity problem at the core of the each DMFT step, then iterated to ap- 
proach self consistency. We terminate the DMFT iteration when Eq. (Q) 
is satisfied for p c to an accuracy of 10 -4 in the relative norm of successive 
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iterates (i.e. ||p™ — /Oc _1 |l/ll/°c II < 10 -4 where \\p\\ = i ^) 2 ); an d the 

chemical-potential self consistency is satisfied for n c to at least 5 decimal 
places, for all fillings and temperatures. 

For the case of a Bethe lattice, used in this paper, the self energy 
So-(u;,T) of Eq. @ can be eliminated from the DMFT self consistency and 
left out of the calculation. However, it is required to calculate quantities 
such as transport coefficients, so we outline here its derivation. We first 
separate out the n = conduction electron orbital and re-diagonalize the 
remaining orbitals, n = 1,2, . . ., to obtain a new set of k-states, ak CT , with 
dispersion E^ a . In terms of these, the effective impurity model becomes, 



+ Y E ^ a i,u a K» + Y (^4, a* + H - c ) ■ ( A - 2 ) 



k,u k,i/ 

The hybridization matrix elements are introduced formally via Ao h\^ a = 
^k a k,o-> which expresses in terms of the new complete set of k-states. 
In order to obtain an expression for the self-energy, we follow the method 

inll, and consider equations of motion for G Cj(T (uj,T) = ((co )Cr ; cj CT )) = 

(w - e Q )G c>a (u,T) = 1 + Y Vw((awy, + T a (u, T), 

k 

(to - £ , k,cr)((ak, -;4,cr)) = V k G c ^(uj,T), 

where 

T a (cj,T) = ((O a ;bl,a)), (A.3) 

and 

O ff = ^(b - <r S-' r + abo^S'). (A.4) 
Solving the above equations for G C)(7 (lu,T) gives 

G c>a (u,T) = [a;-e -A CT (u;,r)-S CT (a;,r)]- 1 (A.5) 

A CT (u;,T) = ( lVk i \ ( A - 6 ) 
^ (w - #k) 



As discussed in Ref. 68, this form of the self-energy is best suited for DMFT (NRG) 
calculations. The "single-particle" renormalization and broadening of the 
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local Wannier orbital due to coupling to the remaining conduction electron 
orbitals n = 1,2,... is contained in the hybridization function A a (u;,T). 
The many-body renormalization and broadening of this orbital is contained 
in the self energy, E CT , which vanishes in the non-interacting limit J — > 0. 
The DMFT self-consistency condition @) reduces in the present case of a 
Bethe-lattice to ^ 

e + A a {u,T) = ^G Ct(T (u;,T). (A.8) 



Appendix B. SPECTRAL FUNCTIONS 

The low-energy part of the /-electron PES of the Kondo lattice model 
(|l]) is determined from the local T-matrix, To, of (|l]), via 

A(u,T) = — 1m.%(u,T), (B.9) 

and Tq is related to the k-dependent T-matrix of the Kondo lattice model 
T k;k /, via 

T (co,T) =^T k>k ,(u;,T) . (B.10) 

k,k' 

7k,k' is defined by the identity 

G k ' k ' = 6 Kk ,G* k + G kk T k)k ,G k ' k '. (B.ll) 

where G kk and G kk are, respectively, the non-interacting and interacting 
conduction electron Green's functions of the Kondo lattice model. An ex- 
plicit expression for Tq has been given in Ref. |69|: Tq = {(O a ;Ot)), where 



the operator O a is defined in Eq. A. 4. Note that A(u, T) only describes 



the part of the /-spectrum lying close to the Fermi level and associated 
with the lattice Kondo resonance. The higher-energy /-electron excitations 
associated with the /-level satellite peaks are projected out in this model. 
Hence, unlike p c , A(cj,T) is not normalized to unity. Its integrated weight 
reflects the weight of the many-body Kondo resonance and will therefore 
depend strongly on temperature. Finally, for the purposes of identifying the 
relevant low-energy scales, it is also useful to consider the local magnetic ex- 
citation spectrum, x" z (u;,T), of the effective quantum impurity model (|3|). 
This is defined by 

X 'L(oj,T) = -- Im Xzz(lo,T), (B.12) 

where Xzzi^jT) is the longitudinal dynamical susceptibility of (|3|) after self- 
consistency is achieved. 
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